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ABSTRACT 

This paper presents a set of new conditions on the augmented density of a spherical anisotropic 
system that is necessary for the underlying two-integral phase-space distribution function to 
be non-negative. In particular, it is shown that the partial derivatives of the Abel transforma- 
tions of the augmented density mu st be non-negative. Ap plied for the separable augmented 
densities, this recovers the result of lvan Hese et alj d201ll) . 
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1 INTRODUCTION 



Recently, Iciotti & Morgantj d20 lOh raised a question whether the 
density slope-anisotropy inequality y > 2p is the necessary con- 
dition for the consistency of the underlying two-integral phase- 
space distribution function. Here, y is the (negative) logarithmic 
density slope whereas fi is the so-called Binney ani s otro py pa- 
rameter. Extending the earlier finding by lAn & Evanj d2006h that 
the in equality is necessary at the centre (given a finite potential 
welD. ICiotti & Morgantii Hoioh have shown that, for wide varieties 
of anisotropic spherical systems built by flexible fami lies of ana- 
lytic two-integral distributions of certain ansatz (e.g., Cuddefordl 



199lllCuddeford & Louis 1995; Baes & van Hesdl2007ir the fore 



mentioned inequality holds everywhere in radial positions given 
that the distribution function is also non-negative everywhere in the 
accessible phase space and that the central anisotropy parameter /? 
is restricted to be fio < \ . They also presented the equivalent con- 
dition to the inequality for the separable augmented density, wh ich 
has subsequently proven bv lVan Hese. Baes. & Deionghe] fefjllh to 
be satisfied by such a syste m if P < ^ . The latter aut hors have also 
shown that the condition o f lCiotti & M organti (2010) is not valid if 
/3o > \, which they have demonstrated with three counterexamples. 

Here, I present new conditions on the augmented density that 
are necessary for the consistency of the distribution function. It is 
found that when they are applied for a separ able augmented density , 
one of the conditions recovers the result of lvan Hese et al.l d201lh . 
providing an alternative proof. 



2 AUGMENTED DENSITY FOR SPHERICAL SYSTEMS 



According to the leans theorem, a general spherical dynamical sys- 
tem in equilibrium is described by the phase-space distribution 
function (DF) of the form of /(£, |L| 2 ) where £ and |L| 2 are the 
two classical isotropic isolating integrals admitted by the spherical 



potential, namely, £ = "F ■ 
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is the specific binding energy and 



L = y \L\ 2 = rv t is the magnitu de of the specific an gular momen- 
tum. Here the notation of lBinnev & Tremainej ( 120081) is followed so 
that "F is the relative potential above the boundary (*F = — G> where 
(t> is the gravitational potential if the system is infinitely extended 
but the potential is escapable) and the bound particles are given by 
£ > 0. In addition, v 2 = v 2 + v 2 where v t and v r are the tangen- 
tial and radial velocities. This paper concerns an escapable system 
or a system with a finite boundary, and therefore the local veloc- 
ity distribution always cuts off (at the escape speed) and £ < is 
inaccessible^ 

Given the DF, its velocity moments are given by 



P„, m (T,r 2 ) = JJj ^dVv 2 "'/(£ = 1 F- y,L 2 = r 2 v 2 ) 



/v 2 "v 2 '" +1 dv t dv r . 



(1) 



Here, the zeroth moment po.oCY, r 2 ) as abivariate function of "F and 
r 2 is referred to as the augmented density. Once the potential ¥(/■) is 
specified, the local density is found by p(r) = pofll^ir), r 2 ] whilst 
all the even-integral velocity moments are given by (v 2 "v 2 "') = 
Pn,mlpa,a with the specified 4* = "F(r). In a self-consistent system 
on the other hand, the Poisson equation with the augmented density 
as the source term reduces to an ordinary differential equation on 
Y, and so p becomes entirely specified by /. 

By changing integration variables, the augmented density 
Po,oC^' r 2 ) is formally expressed to be a bivariate Abel-like integral 
transformation of the DF, that is, 



p ,^,r) 



where 



r 2 IL 



dfidL' 



6>0, L 2 S0,«>0 



f(G,L 2 ) 

% M2 



< K(S,L 2 ;'¥, r 2 ) = 2(Y 



(2) 



(3) 



1 This excludes some systems such as an isothermal sphere. Some of the 
following results are still valid for such systems if they are not dependent 
upon the choice of the lower boundary £ = for £-integral whereas the 
corresponding proper result may be addressed by choosing £ = — oo instead. 
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The integral is over the triangular domain in (£, L 2 ) bounded by 
lines £ = 0, L 2 = and % = 0. The last line is the same as the di- 
agonal line given by £ + L 2 /(2r 2 ) = W. It is known that upon certain 
restrictions on its regularity, equation ((2} can be uniquely inverted 
at least formally to yield the DF, /(£, L 2 ), provided that po.oOP, r 2 ) 
is sufficiently well-be having. Two particular such formulae after 
Hunt er & Oiar] \ 1993b are provided in Appendix iBl 



3 ABEL TRANSFORM OF AUGMENTED DENSITY 

Let us consider two separate changes of integration variables in 
equation l|2}, to (Q,L 2 ) and (S,R 2 ), respectively, where Q = £ + 
L 2 /(2r 2 ) and R 2 = I?/[2(*¥ - £)]. 



3.1 on the potential 

With Q = £ + L 2 /(2r 2 ), it is easy to find that the Jacobian deter- 
minant for the coordinate change (£, L?) — > (Q,L 2 ) is simply the 
unity. Since we then have "K = 2( V P - 0, the augmented density is 
written to be 



f[G m ,AQ\L 2 



popQ¥,r l ) = 2p ff dgdZr 

r2 JJe sm t(0>o,i?x>,(Ki V20F - Q) 



— == dL-/[£ max (0,L-] 
Jo vY-QJo 



(4) 



where £ max (G) = G _ L 2 /(2r 2 ). Here, the inner integral in the last 
line is independent of f, and furthermore, the outer integral is in 
the form of the Abel transform. Consequently, the inner integral 
can be inverted from p 00 by means of the inverse Abel transform. 
That is to say, 



a r 

32 Jo 



e Po,o( > P,r 2 )d4' V2tt 2 r 2r2Q 2 



dL 2 /[£ max (0,L 2 ] 

o 



(5) 



-f 



Next, we can show that 

< 2 Q 



dL 2 f[S m!a (Q),L 2 



rO-P-Q ,2 

2g/(O,2r 2 + J o dL 2 / (1O) [£ max (0,L 2 ] — (6) 



where / (10) (£, L 2 ) = 9//9£. However, we also find that 

a 



32 



dL 2 L 2 /[£ max (0,L 2 ] 



= 2r 



Hence, 



4^/(0, 2r 2 + f dL 2 L 2 / <1-O) [£ max (0, L 2 ] (7) 
Jo 

a r 2,lQ 

- dL 2 /[£ max (0,L 2 ]. 
r Jo 



_9_ 

9T 2 



2 / ° M /7 -4 =^fT dL 2 L 2 /[£ max (0,L 2 ] 
Jo VG- * J V2r 4 Jo 



(8) 

Strictly, the argument so far only indicates that the difference be- 
tween the left- and right-hand sides is independent of Q, but both 
vanish if Q = and thus the equality holds. 



3.2 on the radial distance 

Given R 2 = L 2 /[2(4 / - £)], we first find the Jacobian determinant 
9(£,L 2 ) 



2(f - £) 



9(£,^ 2 ) 

for the coordinate change (£, L 2 ) — > (£, R 2 ) whilst we have that 



(9) 



7C = 2QV - £) • 



2R 2 ( X ¥ - £) 



2QV - £) ■ 



R 2 



(10) 



Then, equation ® reduces to 

2 3 ' 2 7t 



Jo Jo 



d£ , /[fi.L 2 ^ 2 )]. (11) 



Vr 2 -fl 2 " 

Here L 2 lax (i? 2 ) = 2R 2 (*¥ - £), which does not depend on r 2 . The 
integral is now over the rectangular region and so the double inte- 
gral can be performed in any order given the absolutely integrable 
integrand. In addition, the i? 2 -integral is again an Abel transform 
and its inverse transform results in 



_a_ r R 

3* 2 Jo 



rp ,o(y,r 2 )dr 2 
Vfl 2 - r 2 



2 3 ' 2 n 2 f d£V¥Tfi/[£,L 2 lax (^ 2 )] 
Jo 



Next we note that 
r R 2 



Vtf 2 - r 2 



Vfi 2 - r 2 riR 2 -r 2 r 

and therefore (provided that po.odr is integrable over r = 0) 
~ R >A>,o(Y,r 2 )dr 2 



(12) 



(13) 



9 r 

Wo 



Vfi 2 - r 2 
9 \ D 2 C Rl PwW,r 2 )dr 2 
Wi R Jo 



p w QV,r 2 )dr 2 



(14) 



/, p2 3 \ f R >o,o( > I',r 2 )dr 2 
\ 2+R ~w)l Q ,.. 



dR 2 

On the other hand. 



■^R 2 -r 2 



+ 2R 2 \ d6^/¥^Sf (0[) [S,Ll.^R 2 )] (15) 
Jo 

where f°- l \6,L 2 ) = 9//9L 2 assuming lim L ^ L/(fi, L 2 ) = 0. 
However, since 



9 
dR 2 



f[6,Ll^R 2 )] 



equation d 1 5 1 > further reduces to 



2V^£/ (01) [£,L 2 nax (fl 2 )], 



(16) 



9 

9T 



dfiVT^fi/^LL^ 2 )] 

Jo 



( ^ + , 2 ^)r d£ ^p, (17) 



which is in fact valid provided that fdL 2 /L is integrable over L 2 = 
(c.f., AppendixjCjl. Consequently, we finally find that 



_9_ p >o.oCy,r 2 )dr 2 2 3 /2 _ 2 
3 >i 'Jo riR 2 -r 2 



Jo 



(18) 
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4 CONSISTENCY OF DISTRIBUTION FUNCTION 

The results so far are summarized as follows. If we define the Abel 
transformations of the augmented density as 



pC?, r ) = - — = 

n Jo Ry/r 2 -i 



dR 1 



(19) 



pCP,r 2 )_ 

VP 2 ^ 

we find that applying the transform on equation ((2} results in 

V27T 



J J 6, 



d£dZ//(£,Z/); 



E>0, L 2 »0,'A>0 



/(£,£ 2 ) 



(20) 



pCP,r 2 ) = — dfidL 2 £A-J_ 

r JJe»o,i?»o,'K>o l 

pCF, r 2 ) = V2tt ff d£dL 2 

JJsso,l^o,»o VT-£ 

The calculations are performed by interchanging orders of integra- 
tions (so that Q- and A 2 -integrations of the transforms become the 
inner-most ones) and carrying out the Q- and A 2 -integrations of the 
transforms explicitly. The details are given in Appendix lAlfsee also 
IOianl ll993). Not e that the middle equ ation J20t is also identifiable 
to eq. (CIO) of iHunter & Oianl il993l) . 

The partial derivatives of equations \2Q\ are related to sim- 
ple linear integrals of the DF along the diagonal line given by 
£ + L 2 /(2r 2 ) = *P. In particular, 



dpQ¥, r 2 ) V2tt 



3T 



r Jo 



9[r 2 p(T,r 2 )] - '- ; - 

9r 2 VI, 
9pCP,r 2 ) 27t r^ A ,2 



(21) 



>CP,r 2 ) _ 2tt f /r dL , L 2 2 ^ 

„ /"^r 2 * r 2 



9p(f , r 2 ) _ 7t 
9r 2 r- jo 

Here each integral may alternatively expressed through 

~2r 2 ¥ 



-'X 



= (2r 2 )" ,+1 d£ C? - £)"7[£, 2r CP - £)] . (22) 
Jo 

This is verified by changing the dummy integration variables. The 
last line of equations i2U may be expressed in terms of p because 



9pCP,r 2 ) 2 9 v 2 9[rp(T,r 2 )] 



(23) 



g,.2 V2 9r 2 7^ ' ' g,.2 

which may be shown with equation J 14b ■ 

Since the minimal consistency of the DF indicates that f > 
in the accessible region of the phase space and also L 2 > and 
£ < Y, equations d21t are all non-negative, provided that they are 
indeed well-defined. In fact, p and p are well-defined given that 
the DF is integrable - which further implies that po.o( v f> r 2 )dT and 
Po,oCP> r 2 )dr 2 are integrable over *P = and r 2 = 0, respectively. 
However, the definition of p and subsequently the third line of equa- 
tions d21t require /(£, L 2 )dL 2 /L to be integrable over the region 
containing L 2 = 0, which is a strictly stronger condition than the 



mere integrability of the DF. For instance, if /(£, Lr) = ZT 2/? g(£) 
and po.oCP, r 2 ) = r'^AQ¥) with ~ s: B{< 1), it is clear that p cannot 
be defined as in equation U91 whilst the right-hand side of the third 
line of equations J2U diverges. 



4.1 separable augmented density 

Next, we consider the particular cases for which the *P- and r 2 - 
dependences of the augmented density are multiplicatively separa- 
ble, that is, po.oCP^ 2 ) = ACV)B{r 2 ). It is a fairly straightforward 
exercise to establish (c.f, Dejonghe 1986) 



dpi,o 
9Y 

and 

9(r 2 Pi,o) 
dr 2 



Po.o 



f 

Po.o 

Jo 



■ POA 



POA 

2pi.o 



(Q,r 2 )dQ, 



9 log Pi.o 
9 log r 2 



(24) 



(25) 



from equation l[T}. The last result is also directly r elated to the Bin- 
ney anisotropy parameter l lBinnev & MamorJI 19821) . i.e., 



P(r) = 1 - 



<v?) = _ PoAmr),^} = 91og Pl , 
2<v 2 > 2p 1 . [T(r),r 2 ] 91ogr 2 



¥(r),r 2 



(26) 



Therefore, with a separable augmented density, we further have 



Pi,oCP,r 2 ) = S(r 2 ) A(Q)dQ, 
Jo 

which is also separable. In addition, 
a d log B(r 2 ). C2B(r) 



d log r 2 



logS: 



P 



dr, 



which is thus completely specified (up to a constant) given the 
anisotropy parameter. 

The first and last of equations J2U being non-negative for a 
separable augmented density then indicates that 



dACF) 

d*F 
where 



dB(r 2 ) 
dr 2 



•f 

Jo 



A(Q)dQ 



B(r 2 ) = 



f 

Jo 



RB(R 2 )dR 2 
Vr 2 - R 2 ' 



(27) 



(28) 



Here we have also used A(W) > and B(r 2 ) Ss 0, which are direct 
consequences of the DF being non-negative. It is also obvious that 
ACT) > 0, and so we find from the second line of equation j21t that 



d[r 2 B(r 2 )] 
dr 2 



(1-/3)B>0. 



(29) 



Given that B(r 2 ) > 0, this is equivalent to/3 < 1 , which is physically 
obvious. None the less, it is a nontrivial necessary condition on 
B(r 2 ) for the consistency of the corresponding DF. Furthermore, 
we also find that 



dr 2 J 



RB(R 2 ) dR 2 
Vr 2 - R 2 



7-[ 

r Jo 



d[R 2 B(R 2 )] RdR 2 



dR 2 



Vr 2 - R 2 



(30) 



and thus equation ( 1291 ) is in fact a sufficient condition for the second 
part of equation j27\ . Finally, let us suppose that 



Sir 2 ) : 



r 1 B(R 2 )dR 2 _ r 

Jo Rir 2 -R 2 Jo 



S(r 2 sin 2 |)d0 



(3D 
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converges (and so is well-defined). Then since it is also non- 
negative, the non-negativity of the DF indicates that 



dAOF) 



> 0; 



d¥ ' "' dr 2 
Here we note that 



d[rB(r 2 )] _ 1 dB(r 2 ) 
r dr 2 



>0. 



(32) 



(33) 



dAffl A(0) | r dQ dA(Q) 

and so equation J32b is also a sufficient condition for equation 
d27l l. However, equation d32t is the necessary condition for the 
consistency of the DF provided that B{r 2 ) is well-defined in equa- 
tion J3 1 b whereas equation d27t is the necessary condition for the 
non-negativity of any integrable DF that produces a separable aug- 
mented density. 

The partial result of Ivan Hese et al. I {III 1) is also deduced 
from these. If Bo = lim r _,o B(r), this implies that B(r 2 ) ~ r~ 2/S ° 
as r — > 0. Hence, if 2B < 1, then B converges and is well- 
defined. Consequently, equation j32j , in particular, dA/dT > 
is t he necessary condition fo r the consistency of the DF. Follow- 
ing |Qotti&^organt| feOlCt) . this is also equivalent to the density 
slope-anisotropy inequality (y > 2B) for the system with a sep- 
arable augmented density. The result here is however not directly 
applicable for the case that B = \, for which equation J3 1 1 > di- 
verges. 
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APPENDIX A: ABEL TRANSFORMATION OF Pn,„C¥, r 2 ) 

It is easy to establish that the moment function p,,,,,^ 1 ?, r 2 ) defined 
in equation (TJ is all expressed as an integral transformation of the 
DF in (£, L 2 )-space, i.e., 



r lm p,UV, r 2 ) = ^jJ^dSdL 2 <K"-?L 2 '"f(S,L 2 ) 



(Al) 



where the integral is over the domain in the (£, L 2 ) defined to be 
T = j (£, L 2 ) | 8 > 0, L 2 > 0, <K> ). 

Here, L 2 '" = (L 2 )'" is itself a function of only L 2 and so we may set 
m = without any loss of generality in the following discussion - 
formally, F(£, L 2 ) = L 2m f{&, L 2 ), etc. 

Next, we define the Abel transformations of equation dAlt . 



r 2 ) 



If* 

"Jo 



CP- 01/2 ' 



2 i r p„, (v,R 2 )dR 2 

P "^ r ^nJ R^-R^ 2 ' 



(A2) 



p„(T,r 2 ) = - 



1 C" R 2n+l p nfi Q¥,R 2 )dR 2 



Interchanging the orders of integrations indicates that 
Pn (V,r) = 72 JJ T dSdLfJ s ^ W - W2 ; 

[7C(R z )] n ?dR 2 



, rr , r YK(R 2 )f-2dR 2 

P „(>P,0 = 2j|d£dL/j ——j. 



(A3) 



where <K(Q) = K(&,L 2 ;Q,r 2 ) and<K(R 2 ) = 'K{8,L 2 \ X V,R 2 ) with 
% being defined in equation (O. The inner-most integrals can be 
evaluated analytically. In particular, 



r VK(Q)r--dQ <K n 

J 6+ 4 CV-Q) 112 V2 



B(»+±,±); 



i 

X 



YK{R 2 )r-2d^_ _'JC nl , ,, 

1/2 r l V 1 2' 2>> 



_j2_ R?(r 2 - R 2 ) 1 ' 2 rL 

20P-fi) 

r R 2 "- l YK(R 2 )f- l 2dR 2 _ 

{r ^ R 2 } m - V2CP=£) 



(A4) 



B(«+i i), 



where 7C = 7C(£, L 2 ; % r 2 ), and 

1 f-idx r(; 3 +i)r(l/2) (i/2) n7 r 

1) 



B(, 



, i r 1 x"4dx r(n+i)r 

n+2 ' 2 ' " Jo (l-*)" 2 " T(n + 



The last equality is valid for a non-negative integer n, and (a)„ = 
Y\"=o( a + i s me Pochhammer symbol. This may be proven by the 
changes of variables given by 

K(Q) K{R 2 ) _._ 2 . R 2< K(R 2 ) 

x(Q) = -^-\ x(R) = ^-- x(R)= — , (A5) 

respectively. Finally, gathering the results so far results in 

. _ a . (1/2), V2' 



n\ r 2 jj t 



JJ d&dl}% n 



m,L 2 )- 



PnC¥,r 2 ): 



n\ r JJ T L 
d/2)„ 



(A6) 



Y27tr 2 " JJ dfidL^'- i^'f;,^ 



(vp _ 6)1/2 ' 
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APPENDIX B: HUNTER-QIAN INVERSION FOR /(£, L 2 ) 

The Abel transformations of the augmented density examined in 
this paper indicate that the particular problem of inverting the aug- 
mented density for the DF, /(£, L 2 ) is also closely related to finding 
the (even part of) two-integral distribution function / + (£, J 2 ) from 
the axisymmetric density p[ v P(i? 2 ,z 2 ),/? 2 ], as the integral transform 
that relates them can be formally identified with equation J20b . 
Hunte r & Oi ar] jl993h had presented the complete solution for this 
latter problem by means of the complex contour integral, and also 
discussed its relation to the spherical anisotropic case. Although 
they did not provide explicitly the corresponding general inver- 
sion formulae for anisotropic spherical systems, they are in fact 
no m ore complicated than the axisy mmetric cas e faunter & Oianl 
1 19931 eq. 3.1). That is to say (see also loianll 19931) . 

l 9 rV-» dZ -r 



i 9 r 1 

- 2 5/2 7t 2 i g^J £ 



<Tem,)+ 



Z-£ 
dZ 

(Z-£) 1/2 



2(Z-£) 



9\Z, 



(Bl) 



2(Z-£) 



where 



, dpC¥, r 2 ) 



PQ¥, r 2 ) 



9>P 
3p(<F, r 2 ) 



1__9_ f v PofliQ.^dQ 



J__9_ 

7t9T 



r 



CP- 01/2 
PO.qCP,^ 2 )^ 2 

R(r 2 - R 2 ) 1 ' 2 ' 



the below of the real axis, winds along the positive (counterclock- 
wise) orientation, crosses the real axis upwards (from the negative 
to positive imaginary) at the right side of Z = ^^(G), and returns 
to Z = £min from the positive imaginary side. (Note that, given that 
£ < is inaccessible, £ m ;„ = 0.) For more detai l s inclu din g the 
definit ion of Ye^fi), please refer iHunter & Oianl dl993|) or lOiar] 
Jl993h . We note that the formula reduces to the lCuddefordl formula 
if the r 2 -dependence of the augmented density is given by a pure 
power law [ i.e., Pq.qCF, r 2 ) = r _2/J A( v P)]. Of course, this further re- 
duces to the Eddington formula if /? = 0. 



APPENDIX C: ALTERNATIVE DERIVATION OF Eq{l7] 

In equation J 1 5b . we have assumed lim L ^o L/(£, L 2 ) = in order to 
omit the term lirn^g VT* - £/[£, L 2 nax (/? 2 )] in its right-hand side. 
However, this assumption is only incidental here and the result in 
equation (T7J and subsequently equation j 1 8t are valid provided 
that the integrals there are well-defined and differentiable. This may 
be argued based on the fact that equation dl8t is rederived as the 
third line of equations d21t by differentiating equation d20b directly 
and using equation ( 122b . Alternatively, using equation d22t . we first 
find that 



f 

Jo 



d£ -£/[£, 2R 2 (¥ - fi)] 



{2R 2 f 



I 



* /[6,2« a (¥-6)] 1 r 2R ^dL 2 



(2R 2 y 2 



V^-fi 

Then, the direct differentiations indicate that 
9 r* , 

— J o d£V^6/[£,2fl 2 («F-6)] 

i r 2R2 ' 1 ' 1 2 

= >F 1/2 /(0, 2R 2 V) + — — di?Lf lfi Hy - At, l 2 \ 

J (2R 2 ) 3 I 2 J J V 2R 2 /' 



and 



and pCP, r 2 ) and pOP ', r 2 ) are as defined in equation d 1 9b . Here, the 
outer integrals of equation JB 1 b are contour integrals in the complex ^ 
Z-plane whose path is given such that it starts from Z = fimin to 



. Jo 



¥ f[S,2R 2 W-S)] 
dfi- 



= /(0, 2R 2 V) + — — - dL 2 — /< 10 )(T - 

R J 2 3 ' 2 R 4 J L V 

That is to say, 



9 



Jo 



T /[£,2tf 2 0F-£)] 
d£- 



VT-fi 

1 g pf 

= -g^J o dfiV^£/[£,2i? 2 CP-£)], 
which is equivalent to equation ( 117b . 
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